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Radical-ion-pair reactions were recently shown to represent a rich biophysical laboratory for the 
application of quantum measurement theory methods and concepts. We here propose a concrete 
experimental test that can clearly discriminate among the fundamental master equations currently 
attempting to describe the quantum dynamics of these reactions. The proposed measurement based 
on photon statistics of fluorescing radical pairs is shown to be model-independent and capable of 
elucidating the singlet-triplet decoherence inherent in the radical-ion-pair recombination process. 



I. INTRODUCTION 

Spin-selective radical-ion-pair reactions represent a 
rich biophysical/biochemical system in which spin de- 
grees of freedom can dramatically influence the fate of 
biologically significant chemical reactions. The study of 
radical-ion-pair reactions is at the core of spin chemistry 
PI, by now a mature research field directly related to 
photochemistry Q and photosynthesis Q. Radical-ion- 
pair reactions determine the late-stage dynamics in pho- 
tosynthetic reaction centers 4, 5], and furthermore there 
is increasing evidence that radical-ion-pair reactions un- 
derlie the avian compass mechanism, i.e. the biochemical 
compass used by migratory birds to navigate through the 
geomagnetic field |6l— ll 21] . Hence the fundamental under- 
standing of these reactions if of high scientific interest. 

Radical-ion pairs are molecular ions created by a 
charge transfer from a photo-excited D*A donor-acceptor 
molecular dyad DA, schematically described by the reac- 
tion DA -t D*A -> D* + A*~, where the two dots rep- 
resent the two unpaired electrons. The magnetic nuclei 
of the donor and acceptor molecules couple to the two 
electrons via the hyperfine interaction, leading to singlet- 
triplet mixing, i.e. a coherent oscillation of the spin state 
of the electrons. The reaction is terminated by the re- 
verse charge transfer, resulting to the charge recombi- 
nation of the radical-ion-pair and the formation of the 
neutral reaction products. It is angular momentum con- 
servation at this step that empowers the molecule's spin 
degrees of freedom to determine the reaction's fate: only 
singlet state radical-ion pairs can recombine to reform the 
neutral DA molecules, whereas triplet radical-ion pairs 
recombine to a different metastable triplet neutral prod- 
uct. 

The fundamental quantum dynamics of radical-ion- 
pair (RP) reactions rests on a master equation satisfied 
by p, the density matrix describing the spin state of the 
molecule's two electrons and magnetic nuclei. This mas- 
ter equation has to describe (i) the unitary evolution of 
p due to the magnetic interactions within the radical-ion 
pairs, which is straightforward, (ii) the loss of radical-ion 
pairs due to the recombination reaction leading to the 
creation of neutral products and (iii) the state change 
of unrccombined radical-ion pairs. The perplexity of the 
combined presence of all those phenomena is partly the 



reason behind the ongoing debate on the particular form 
of this master equation. The current standing of this 
debate is the following. Kominis derived 13[ a mas- 
ter equation for the term (iii) and put forward a master 
equation [l4j taking into account the reaction term (ii). 
Another master equation was put forward by Jones and 
Hore [l5j . while several authors fl6l - [l8j argued in favour 
of the traditional master equation of spin chemistry. So 
the same physical system is currently described by three 
theories. This situation is clearly unsatisfactory, and al- 
though theoretical arguments could in principle point to 
the fundamentally correct theory [19|, l20(, the need for 
an experiment with discriminatory power is obvious. We 
will here propose exactly such an experiment. 



II. PHOTON STATISTICS 

We will here describe the proposed measurement and 
explain the physics behind it in the next section. We con- 
sider radical pairs (RPs) that recombine only through 
one channel e.g. the singlet, and hence &t = 0. We 
also consider each recombination event to be accom- 
panied by a photon emission through e.g. an exci- 
plex fluorescing molecule. Let n t and n t +dt represent- 
ing the photon counts in the time intervals (t, t + dt) 
and (t + dt,t + 2dt). The stochastic variables n t and 
rit+dt have expectation values Nt = ksdtTr{Qsp t } and 
N t +dt = ksdtTr{Q s p t+dt }, respectively, given by the sin- 
glet recombination products during the respective time 
intervals. The probability to actually observe nt {nt+dt) 
photons is given by the Poisson distribution with expec- 
tation value Nt (N t +dt)- 

We now define the stochastic variable Sn = (n t +dt — 
n t )/N t . If n\ and n% are Poisson stochastic vari- 
ables with mean values N\ and N%, respectively, then 
the probability that the difference n\ — n<i = k fol- 
lows the Skellam distribution, given by f(k;Ni,N 2 ) — 
e-^ Nl+N2 \N 1 /N 2 ) k/2 I\ k \(^VW^2), where I k (x) is the 
modified Bessel function of the first kind. The mean and 
variance of the Skellam distribution are p = N\ — N 2 and 
a 2 = Ni + N 2 . Therefore the mean and variance of Sn is 
Sn = (Nt+dt - N t )/N t and of„ = (N t + N t + d t)/N?, re- 
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FIG. 1: We numerically integrate the Jones-Hore and the 
Kominis equation for a magnetic Hamiltonian of the form 
H — cjiSiz + uj2S2z- Starting with 10 12 RPs, we then calcu- 
late the expected photon counts in a time interval At from 
N * = It +At ks{Qs)dt, and form the ration (N t - N t+A t)/N t , 
which is plotted for the two theories. We also plot the sta- 
tistical error a, almost the same for both theories. We chose 
At = 1/Aks in order to have enough statistics for a to be 
small enough to allow a statistically meaningful comparison 
between data and theoretical expectation. 



arbitrary hyperfine interactions. At high enough mag- 
netic fields, we ignore all hyperfine interactions and keep 
just the Ag contribution to S-T mixing. Thus the theo- 
retical comparison between the two theories is completely 
robust and independent of the details of the molecule's 
hyperfine interactions. In Fig. la we plot the expectation 
value of Sn according to the theory of Kominis (Snk) and 
Jones-Hore (Sujh)- We also plot the standard deviation 
asn (roughly the same for both theories) . It is clear that 
the oscillations in Snx decay rather slowly in time. In 
contrast, the oscillations in Sujh decay fast, with a time 
constant on the order of the recombination rate k$- Fur- 
thermore, the statistical error ag n remains small enough 
for the oscillations in Sn/, to be detectable for a long time 
well beyond l/ks- 



III. PHYSICAL EXPLANATION 

The expected number of photons detected in the time 
interval between t and t + dt will be N t = kg(Qs) t dt, 
whereas those detected in the following time interval be- 
tween t + dt and t + 2dt will be N t +dt = ks(Qs)t+dtdt. 
Hence Sn = d{Qs)t/{Qs)t, where d(Q s )t = (Qs)t+dt - 
(Qs)t- Obviously, Sn is maximum when the slope 
d(Qs)t/dt is maximum. The maxima of the slope 
d(Qs)t/dt occur in between the peaks and troughs of 
(Qs)t, which are points of minimum singlet-triplet (S-T) 
coherence. In other words, the maxima of Sn occur at in- 
stants in time when there is maximum S-T coherence. As 
shown in Fif.lb, where we depict the expectation value 
(Qs)t derived from the two theories, the maxima of the 
slope d(Qs) t /dt stay roughly constant for Kominis' the- 
ory, whereas they quickly decay for the Jones-Hore the- 
ory. This slope represents the "amount" of S-T coher- 
ence, and hence it is this point where the two theories 
fundamentally differ. We will now analyze the root of 
this difference. 



A. Singlet-Triplet Coherence 

Multiplying p from left and right by 1 = Qs + Qt, 
it follows that p can be written as p = p + p, where 
P = QspQs + QtpQt and p = QspQt + QrpQs- The 
incoherent part of p is represented by p, whereas the S-T 
coherent part is p. From the master equation for dp/dt 
the rate of change dp/dt can be readily calculated. 

B. Prediction of the Jones-Hore Theory 

The Jones-Hore master equation is (for the considered 
case of having kr — 0) 

dp/dt = -i[H,p]-k s (p-Q T pQ T ) (1) 

leading to 

dp/dt = —ksp~ + unitary terms, (2) 

where "unitary terms" =— iQs[H, p]Qr — iQr\H, p]Qs- 
So the coherence p decays away at a rate ks- 

C. Prediction of Kominis Theory 

The master equation derived by Kominis reads for this 
case (fc T = 0) Qj 

dp/dt = -i[H )P ]-JL (q sP + pQ s - 2Q sP Qs) 
- (1 - Pcoh)ksQspQs 

dn s p , . 

Pcoh dt Tr{p} {6) 
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where dns — ksdt(Qs)t andp co h is the measure of 
S-T coherence introduced in [14]]. The second, trace- 
preserving Lindblad term of (j3|) takes into account the 
S-T decoherence brought about by the continuous mea- 
surement of the RP's spin state induced by the intra- 
molecule reservoirs [l3|, [3] ■ The two last terms are the 
reaction terms that change p due to the RP's that recom- 
bined into neutral products in the time interval (t, t + dt). 
This master equation leads to 
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unitary terms (4) 



Neglecting for the moment the second term in the paren- 
thesis, it is seen that the S-T coherence fades away at 
half the rate in this theory compared with the Jones- 
Hore theory. However, this is not the only reason for the 
faster decay of the amplitude of the S-T oscillations in 
the Jones-Hore theory evidenced in Fig. lb. To explain 
why, we rewrite the Jones-Hore master equation as 

dp/dt = -i[H,p\- ks(QsP + pQs - ZQspQs^J 

- ksQspQs (5) 

Comparing the respective terms in ([5]) and ^ , it is seen 
that (i) as noted before, the S-T decoherence takes place 
at double the rate in the Jones-Hore theory, and (ii) set- 
ting p co h = for all times in (|3]) we retrieve the reaction 
term of the Jones-Hore theory. So the amount of spin 
coherence, embodied in the relative balance of p versus 
p (and quantified by p C oh), is not only diminished faster 
due to the first term in ([5]), but also due to the fact that 
the reaction term in ([5]) distorts this balance. This point 
is related to the way the two theories deal with the up- 
date of the density matrix p due to recombination, when 
p is a state having large S-T coherence. This point has 
been analyzed in some detail in [3], but we will here 
reiterate this analysis in the present context. 

At high magnetic fields the T± states are out of reso- 
nance with the singlet, so we can assume that S-T mix- 
ing occurs just between \S) and \Tq). Ignoring hyper- 
fine couplings, we effectively deal with a two-level sys- 
tem. Suppose that at time t the spin state of the RP 
is the pure state \ip t ) = \S). This state has obviously 
zero S-T coherence. Suppose also that after the elapse 
of time r the state \S) is transformed by the Ag mech- 
anism to \ip t +r) = {\S)\ + \T ))/V2. The correspond- 
ing density matrices will be p t = |>5)(S , | and pt+ T = 
(\S){S\ + \T Q ) (T 1 + \S) (T | + \T ) {S\)/2. Let us neglect for 
this discussion the change of the density matrix of the un- 
recombined RPs, and focus on the change of the density 
matrix due to the recombining RPs. The update of pt due 
to RPs recombining between t and t + dt is the same for 
both theories, since this is the incoherent limit. However, 
the update of pt+ T due to RPs recombining between t + r 
and t + t + dt differs fundamentally in the two theories. 
In the Jones-Hore theory pt+ T changes due to recombi- 
nation by dp r — —ksdtQsPt+rQs — —ksdt\S)(S\. In 
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FIG. 2: Figure 4 of [21J. Photon counts of singlet recombining 
(diphenylsulphide-dio) + /(p-terphenyl-di4)~ at various mag- 
netic fields, increasing from 1 to 4. At the highest magnetic 
field one observes quantum beats stemming from the Ag mix- 
ing mechanism. The "calibration" trace 1 at the lowest mag- 
netic field is seen to be governed by two rates, a fast rate on 
the order of ks and a slow rate stemming from the strong 
projection produced by ks on the state that has oscillated to 
the triplet. In any case, the quantum beats decay at a rate 
much smaller than ks, which is what the theory of Kominis 
predicts. 



contrast, in the Kominis theory dp r = —dnspt+r, where 
dns — ksdt(Qs)t+r is the number of RPs that recom- 
bined within the time interval (t + r,t + r + dt). Stated 
qualitatively, in the theory of Kominis and in the case of 
maximal S-T coherence, one removes the complete den- 
sity matrix pt+r as many times as many RPs recombined. 
In contrast, the Jones-Hore theory projects out of pt+ T 
the singlet part QsPt+rQs, no matter how coherent the 
state pt+T is. This has the result that the rest, surviving 
radical-ion pairs appear to be less coherent. In contrast, 
the surviving RPs in the Kominis approach have not lost 
any coherence at all. The artificial loss of S-T coherence 
due to the projective recombination term in the Jones- 
Hore theory is thus the main reason behind the rapid 
decay of the oscillations seen in Fig. lb. 



IV. SPIN RELAXATION 

The obvious difference between the predictions of the 
two theories will fade away if other spin-relaxation mech- 
anisms arc dominant beyond the S-T decoherence process 
internal in the molecule. We emphasize that the latter 
decoherence mechanism, present in both theories, is an 
unavoidable one fundamentally linked to the charge re- 
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combination process. It 's p hysical origin has been ex- 
haustively explained in [l3j and 14]. Other relaxation 
mechanisms, for example spin-changing collisions with 
other molecules, can in principle be suppressed, e.g. by 
performing the measurement at low enough tempera- 
tures. In any case, if the spin relaxation rate of all such 
other mechanisms is denoted by k sr , we can state that 
that in order for the difference between the two theo- 
ries to surface, k sr must be significantly shorter than ks- 
Otherwise, S-T coherence will be quickly suppressed, i.e. 
Pcoh — > 0, and then the two master equations will yield 
qualitatively similar results, as the respective reaction 
terms will coincide. The condition k sr < ks can be ex- 
perimentally satisfied. Indeed, there already exist exper- 
imental data that provide evidence for our assertions. In 
21] Molin presents photon counts of singlet recombin- 
ing (diphenylsulphide-dio) + /(p-terphenyl-di4) _ at vari- 
ous magnetic fields, demonstrating the oscillations due to 
Ag mixing at high fields. The data is reproduced in Fig. 
2. It is evident that there exist two rates a fast rate (on 
the order of ks) and a slow rate (much smaller than ks)- 
It is also evident that while the RP population decay is 
governed by the fast rate, the quantum beats decay at 



the slow rate. This is in contrast to the predictions of 
the Jones-Hore theory, according to which the amplitude 
of the quantum beats should decay at the rate ks- 



V. CONCLUSIONS 

It is important to stress that the particular measure- 
ment we propose is model-independent, in the following 
sense. Clearly the time evolution of observables like (Qs) 
or the magnetic-field effect (MFE) are predicted to be dif- 
ferent by the two theories. However, in practice it would 
be almost impossible for the measurement to discern ab- 
solute differences in (Qs) or the MFE signal, since it 
would be possible to attribute those to an imperfect un- 
derstanding of e.g. the RP's magnetic interactions. In 
contrast, in the measurement we propose, normalizing 
the photon difference by N t largely alleviates this prob- 
lem, and the two theories predict a clearly different trend, 
no matter what the details of the magnetic interactions 
are, i.e. how many nuclear spins there are or what are 
the exact values of the hyperfine couplings. 
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